Shock waves and Birkhoff 's theorem in Lovelock gravity 



E. GravaniiQ 
(Dated: September 13, 2010) 

Spherically symmetric shock waves are shown to exist in Lovelock gravity. They amount to a 
change of branch of the spherically symmetric solutions across a null hypersurface. The implications 
of their existence for the status of Birkhoff 's theorem in the theory is discussed. 
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I. INTRODUCTION 

In general relativity Birkhoff 's theorem is roughly the 
statement that outside a spherically symmetric (even 
time- varying) source the metric field is necessarily static 
and it is given by the Schwarzschild solution [l,][2][3]. 
A better statement of the theorem is that the vac- 
uum field equations imply that a spherically symmet- 
ric C 2 solution is locally equivalent to a maximally ex- 
tended Schwarzschild metric H 1 . An important refine- 
ment came with the work of Ref. Q, building on the 
work of Refs. [8] [9]: Schwarzschild metric is the unique 
spherically symmetric family of vacuum solutions and 
Birkhoff 's theorem holds even if we lower differentiabil- 
ity class to C°. The result which encapsulates, one may 
say, the essence of what makes the theorem possible, is 
that spherically symmetric shock waves do not exist. In 
more words, there cannot be a vacuum solution respect- 
ing spherical symmetry whose first derivative becomes 
discontinuous across a null hypersurface. This result 
was obtained already in [9] . Had such a solution existed, 
any statement of Birkhoff 's theorem for smooth metrics 
would be strongly weakened. On the other hand, the in- 
existence of the spherically symmetric shock waves may 
be attributed to the uniqueness of the spherically sym- 
metric smooth metric of a given mass. There simply 
cannot be a non-trivial matching across a (null) hyper- 
surface respecting the symmetry. 

In dimensions higher than four, Einstein field equa- 
tions can be formally thought of as a special case of the 
Lovelock field equations [lfj. The additional terms are 
of higher order in curvature nonetheless they still arc 
second order differential equations for the metric ten- 
sor. Originally, at least in the recent years, an inter- 
est in Lovelock gravity came through studies on low en- 
ergy effective actions in string theory [ll|. In the last 
decade the interest was revived due the higher dimen- 
sions becoming popular, essentially through works such 
as [lH [l3| . In the last few years a certain amount of in- 
terest in Lovelock gravity was further drawn, taking also 
an attractive turn. The shear viscosity/entropy density 
bounds for a conformal field theory living on the bound- 
ary of AdS 14|[l5| are modified and restricted when ap- 



plying causality conditions on the conformal field theory 
with a Lovelock gravity dual in the bulk [l6|[l7|. More- 
over these restrictions coincide unexpectedly to ones ob- 
tained through positive energy conditions on scattering 
processes in super-conformal field theories 
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findings led to a series of new works, see e.g 
related to Lovelock gravity. 

Establishing Birkhoff's theorem in Lovelock gravity 
requires a more careful phrasing than in general rela- 
tivity. Perhaps the most complete presentation of the 
theorem has been given in Ref. [3l| building on the work 
of [HJ . Comments in that direction have been presented 
in other works [33l - l36l ]. The statement of the theorem 
is one of uniqueness: The vacuum spherically symmetric 
metrics of differentiability class C 2 and for generic val- 
ues of the couplings of the theory, are locally equivalent 
to a specific family of solutions, some of which are the 
Lovelock black hole metrics. A characteristic feature of 
Lovelock gravity is that these solutions are multi- valued. 

Now, what about metrics of differentiability lower 
than C 2 ? In Ref. [37j spherically symmetric vacuum so- 
lutions with C° piecewise C°° metrics where explicitly 
constructed. Their metric is everywhere smooth except 
at certain time-like as well as space-like hypersurfaces 
where it is only continuous; its first derivative is discon- 
tinuous. That is, non-trivial vacuum shells with time- 
and space-like trajectories were explicitly shown to ex- 
ist. These imply two things. One, there is a problem 
of non-causal evolution due to the spacelike vacuum dis- 
continuities. Second, any staticity interpretation of a 
Birkhoff's theorem in this theory is very much weakened 
by the possibility of a series of spherical vacuum shells 
whizzing along spacetime. The existence of spherical 
shock waves in Lovelock gravity, that is of spherical vac- 
uum shells moving along light-like trajectories, are the 
subject of the present work. They elegantly appear as a 
mere change of branch along a null hypersurface. Now 
even if one finds a way, by some nice principle, to censor 
out the non-null discontinuities, the null ones cannot be 
excluded. They are gravitational shock waves, perfectly 
existent in general relativity, with the exception of the 
spherical ones which gives rise to Birkhoff's theorem. In 
Lovelock gravity the spherical ones exist and appear in 
fact as fairly natural objects in the theory along with its 
multiple branch solutions. Presumably, all these non- 
trivial vacuum shell spherical symmetric solutions can 
be regarded as large (non-linear) scalar perturbations of 
a given smooth spherically symmetric metric. 

In order to show that configurations such as the ones 
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discussed in the previous paragraph are indeed solutions 
of a given theory, requires to formulate and solve, one 
way or another, junction or matching conditions in that 
theory. In general relativity, the matching conditions for 
the time- or space-like hypersurface where given their 
final form in Ref. [H]]. Shock waves and null hypersur- 
face matching conditions have been studied in various 
old and more recent works (8| (39l - l47j . The problem we 
set ourselves to solve is to construct certain shock wave 
solutions in Lovelock gravity. To do this we shall ex- 
ploit the formulation of matching conditions first pre- 
sented in Ref. 48] and further elaborated in [49( and 
[Hcj . The basic facts of the method can be summarized 
as follows. The gravitational field is described by the 
vielbein and the spin connection. Consider a sequence of 
smooth configurations which are arbitrarily close approx- 
imations of a given discontinuous one, e.g. one where the 
spin connection becomes discontinuous across hypersur- 
faces [which may also intersect]. The hypersurfaces di- 
vide spacetime up into 'bulk' regions where fields are 
smooth. At the hypersurfaces the discontinuous fields 
are ill-defined. Then the action functional of the theory 
evaluated on that sequence can be shown to converge to 
a new action functional involving only well defined in- 
formation: fields outside the hypersurfaces. The fields 
in the vicinity of each side of the hypersurface contain 
all the information about the discontinuity, therefore 
nothing more needs to appear in the equations. [That 
presumably means that no explicit reference to the in- 
duced fields on the hypersurface is required.] Working 
this way the hypersurfaces need only be locally smooth; 
their causal character, null or non-null, is irrelevant. The 
whole analysis is done off-shell, therefore the equations 
of motion for the discontinuous fields obtained by Eulcr- 
Lagrange variation of the new action are the well defined 
limit of the usual equations of motion for the smooth 
fields. The new action contains explicitly terms for each 
hypersurface. The matching conditions are the parts 
of the equations of motion with support at the hyper- 
surfaces. The convergence of the off-shell action to the 
new action makes them well defined: the bulk fields on 
the sides of a hypersurface of any smooth configuration 
approximating well the discontinuity must obey those 
matching conditions. 

The way of 'matching' in the formulation of 48~5(J is 
somewhat different than usual: One works with globally 
defined fields with a certain amount of discontinuity at 
some places and attempts to write down an action func- 
tional for these fields which is well defined in the sense 
stated above. The spirit of the analysis reflects nicely the 
old work of Papapetrou and Treded Q in general rela- 
tivity we already mentioned. The principle is that field 
configurations of differentiability lower than the order 
of the field equations, are also solutions of the theory if 
they can be approximated arbitrarily well by sequences 
of smooth solutions. A subtlety here is that such se- 
quences is not always easy to make explicit. One such 
case is the vacuum solutions we are interested in. On the 



other hand the well-defined-ness of the matching condi- 
tion is established by sequences of smooth configurations 
off-shell. There are always non-vacuum smooth config- 
urations approximating arbitrarily well a given vacuum 
discontinuous configuration. In that limiting sense so- 
lutions of the matching conditions may be regarded as 
meaningful stationary point of the classical action. 

As the formulation of the matching conditions we shall 
use is not a well known one, we will devote the section [TTT] 
in explaining the technical and physical details leading 
to the matching condition, equation (IT51) . In section HVl 
the formulation is suitably applied to obtain the alleged 
shock wave solution, and certain additional comments 
are given in the final section. 

II. LOVELOCK GRAVITY 

To fix ideas consider a specific example of Lovelock 
gravity. Let our theory be Einstein gravity with cosmo- 
logical constant A supplemented by the quadratic Love- 
lock term in five dimensions. This is usually termed as 
Einstein-Gauss-Bonnet gravity. We denote the coupling 
of the quadratic Lovelock term by a. It has length di- 
mension L +2 . The Lagrangian of the theory reads 

S=-\ [ ^—g+ (I) 

J M 

/ V=g{R + a(R 2 - ±R aP R afi + R a p l6 R a ^ 5 )} . 
Zk Jm 

The spherically symmetric C 2 solutions of the theory 
(JTJ are the Boulware-Deser metrics [5ll - l54j (their general 
forms in Lovelock gravity where first studied in |55|): 

ds 2 = -fBDdt 2 + j B x D dr 2 + rW , (2) 
r 2 [ I 4an 2 A 16aM] 

^ = 1 + s( 1± V 1 + — + — }• 

dil 2 is a metric of the unit round 3-sphere. M is a con- 
stant of integration related to the mass of the solution. 
The Boulware-Deser metrics are not single-valued: they 
have two branches, corresponding to the ± sign choice, 
which have very different properties. For one thing, there 
are no black hole solutions in the metrics of the branch 
corresponding to the sign +. Also for A = this branch 
of the solution is not asymptotically flat: it is asymptot- 
ically anti-de Sitter for a > and de Sitter for a < 0. 
A related fact is that only the metrics of the branch — 
converge to Einstein gravity metrics in the limit a — > 0. 
It has been argued that the asymptotic vacuum of the 
branch + metrics is unstable due the ghost excitations 
[f5lj][56j]. On this basis one may discard it from the be- 
ginning, but that is rather too quick. These metrics are 
an inherent property of Lovelock gravity and their ef- 
fects should be studied before they are discarded with 
reason, or not discarded at all. For one thing, their 
presence might make things better or worse, stability- 
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or other-wise, for the well behaved branch metrics; both 
situations are of physical importance. One should bear 
in mind that the two types of metrics are naturally in- 
termingled by the theory. Lovelock gravity allows C° 
piecewise smooth non-trivial vacuum configurations to 
be cut-and-paste constructed out of metrics of the same 
or different branch as shown in Ref. [37|, see also 57 1. 
New such configurations associated with shock waves are 
presented here. 



III. THE LAGRANGIAN AND EQUATIONS OF 
MOTION 

A. First order formalism 

Proceeding with the mathematical analysis, it is far 
more convenient than working with the metric tensor 
g^u directly to use the vielbein E ° and an 50(4, 1) 
spin connection ui^ b one-forms over M as our variables. 
They will be treated as differential forms E a = dx^E^, 
uj ab = dx^uj.° ,b . We will do our diffentiations and in- 
tegrations using exterior calculus 2 . The metric tensor 
Vab = ( — h • • • +) and the volume anti-symmetric tensor 
(form) e a bcde are the invariant tensors 50(4, 1). From 
them we can build invariant forms over the spacetime 
manifold. 

We will use a convenient notation for contraction with 
the volume form, for example 

e(flE 3 ) := e abcde n ab AE c AE d AE e . 

Its convenience can be seen in variations and differenti- 
ations, e.g. 5{e{nE 3 )} = e(5flE 3 ) + 3e(flE 2 SE). The 
invariance of the contraction also implies that 

d( e (...))=e(D(...)), 

where D is the covariant derivative associated with the 
spin connection. From now on the wedge symbol A will 
be dropped in the wedge product of forms. 

The spin connection throughout this work will be Levi- 
Civita i.e. torsion-free: T a := DE a = dE a +ui a b E b = 0. 
d is the exterior calculus derivative operator. Its nilpo- 
tence, dd = 0, holds on all C 2 forms. Then the affine 
connection on the tangent bundle is Levi-Civita and de- 
fines the usual covariant derivative of general relativ- 
ity. The spin connection is not a tensor. The curvature 
il a b := duj a b +uj a c uj c b is an 50(4, 1) tensor and a two-form 
over M. We can write that more compactly treating the 
forms ui and £1 as matrices: f2 = du> + lo 2 . The curvature 
satisfies the Bianchi identity: DO, = dil + ujfl — ftui = 
identically, as one may verify. 

The curvature form is related to Riemann tensor by 



the relation n ab = \E*E b W v p(J dx p dx a . In other 



words to obtain the Riemann tensor form the curvature 
form one has to invert the matrix Ef. We shall not 
need to do that as we will work exclusively with forms, 
though invcrtibility will not be evaded anywhere. Thus 
the veilbein and connection formulations is equivalent 
to the metric formulations. Then the action (JlJ can be 
written in the following form 

S=-\c a [ e(£ 5 ) + ici / e(f!£ 3 ) + i C2 f e(fifiS) . 
1 Jm 1 Jm 1 Jm 

(3) 

The coupling constants c , c\ and c 2 are related to the 
more usual couplings by 
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a 

72 



(4) 



and k 2 — 8ttG, where we introduce G as the Newton's 
constant (though conventions differ in the literature). 
The cosmological constant appears usually as A = k 2 X. 
We shall prefer using mostly the c couplings. The length 
dimension of c± is L~ 3 and the dimension of C2 is L~ x . 



B. Manufacturing a discontinuity 

Imagine then spacetime consisting of three regions: 
two open bulk regions N, N and a closed region A sur- 
rounding the hypersurface E we would like to be locus 
of the discontinuity. The action may written as 



5= / £+ f C+ f C, 
Jn Ja Jn 



2 Convenient references for this formulation arc 581 and [59 



where £ is a Lovelock gravity Lagrangian. 

Consider the limit A — >• E. It is not that this part of 
spacetime actually shrinks but rather the values of the 
fields E and uj are moved. This can be done by choosing 
a sequence of configurations imitating this process. On 
the common boundary of A with each of N and N the 
values of the fields E and uj are held fixed (and in general 
different). These values become the values of the fields 
on each side of E as it is embedded in each of TV and 
N. These are bulk fields and they are well defined as 
long as there is a coordinate neighborhood covering the 
hypersurface E. 

The limit we described here as A — > E, or any limit 
equivalent to it, was studied for any theory built out of 
forms fields alone, such as E and uj, in Ref. 50]. The 
general result is as follows. Evaluate the action 5 on an 
arbitrary sequence of configurations consistent with the 
limit A — > E as described above. Then for each such se- 
quence of configurations one obtains a sequence of func- 
tionals, whose convergence needs to be checked. Conver- 
gence depends on the amount of discontinuity allowed at 
E. That should be formalized into specific conditions on 
the fields, the 'continuity conditions'. Choosing them ap- 
propriately the action converges to a new action, which 
necessarily depends only on well defined quantities, the 
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bulk fields. [That convergence we usually describe as 
'well-defined-ness' of the new action.] Such continuity 
conditions are not unique and depend on the theory; in 
general they can be quite strange. In the limit A — > E 
the action necessarily involves distributional fields; on 
the other hand the new action depends only on the fixed 
bulk fields on the sides of E. Everything is done off-shell, 
therefore whatever holds for the actions it holds also for 
the equations of motion. The convergence of the ac- 
tion to the new action means then two things. First, the 
equations of motion of the distributional fields are mean- 
ingful, and secondly, they have been neatly disentangled, 
the necessary integrations having been essentially done, 
into the form of the equations of motion deriving from 
the new action. One may then say that the new action 
is an equivalent form of the action of the theory for dis- 
continuous fields. 

A straightforward application of these ideas is offered 
by Lovelock gravity under the most usual continuity con- 
dition. In 50] it was shown that that one may write 
down a new action for Lovelock gravity in the presence 
of discontinuities if the vielbein E is continuous and the 
spin connection co is discontinuous, a result essentially 
contained also in . Lorentz transformations are still a 
local symmetry of the new action [48| [49( . Therefore our 
continuity condition explicitly is: The vielbein is contin- 
uous across the hypersurface modulo a Lorentz transfor- 
mation. This continuity condition will be applied and 
explained in more detail in section IIVI 

Based on these results we may write the action ^ in 
the presence of the discontinuity at E in the equivalent 
form: 

S = -\c, J e{E 5 ) (5) 
+ i Cl J e (^ 3 ) + i Cl J^((u-uj)E 3 ) 
+ ic 2 J e(QQE) 

+ ic 2 J c((w - w){(fl + n) - ~(w - uj) 2 }E) . 

The bulk integrals are over N[jN. Writing (Q — u) 2 
matrix multiplication of the forms is understood. u> and 
f2 are respectively the spin connection and its curvature 
on ./V evaluated at E; more precisely, as N is open, they 
are the limits of these variables as we approach E from 
within AT. Similarly for u> and H in N. The veilbein is 
continuous modulo a Lorentz transformation therefore 
we can express everything in terms of a single veilbein 
field E. 



C. Equations of motion 

Equations of motion are obtained by varying the ac- 
tion with respect to the fields, E and to. 



It was shown in [48[ and more elegantly in [49| that 
the variation of the action w.r.t. the spin connection 
u> trivially vanishes upon imposing that torsion is zero 
i.e. that the connection provides the familiar Levi-Civita 
covariant derivative and ([3]) is indeed equivalent to the 
Lovelock gravity action ([T]). That is, all equations of mo- 
tion, including the matching conditions, are obtained by 
varying w.r.t. the vielbein E. This is a merely algebraic 
variation of the Lagrangian. 

Explicitly, equations of motion are obtained for every 
smooth submanifold in the problem, that is, the bulk 
regions, the hypersurfaces, and in the general case the 
intersections of the hypersurfaces. There is a term in 
the Lagrangian for each one of these submanifolds. Let 
5E a = X(x) a b E b be the variation of the vielbein, where 
X(x) a b is an arbitrary smooth field with support in any 
one of the submanifolds. The equations of motion for 
that submanifold are obtained from the corresponding 
term in the Lagrangian. From the bulk regions we ob- 
tain the usual smooth field equations of motion; from the 
hypersurfaces and their intersection we obtain a set of 
matching conditions. If in general there are other fields 
with energy tensor T a b , the field equations of each sub- 
manifold read 

SC = T a b X(x) b a (6) 

times its volume element. 
In vacuum, T a b = 0. That is, 

SC = . (7) 

In vacuum there is no need to identify components; one 
simply sets to zero all non-trivial terms in this equation. 

Here, we have the bulk regions N, N and the hyper- 
surface E to deal with. 

The field equations in the region N read 

- 5c e(E i SE) + 3ci e(QE 2 SE) + c 2 e(W5E) = . (8) 

A similar expression holds in N for the barred fields. 
Consider the spherically symmetric metric 

ds 2 = -g 2 dt 2 + ^ + r 2 dil 2 , (9) 

where g = g(r) and dil 2 is the metric of the unit round 
3-sphere. The metric (j9]) can be written in terms of the 
vielbein one-forms 

E°=gdt, E x = —, E i = rE i . (10) 
9 

E % is a veilbein of the unit round 3-sphere. We denote 
by ui tJ its Levi-Civita connection. The curvature of a?- 7 
necessarily is fi 4 - 7 = E l EK The Levi-Civita connection 
uj ab of E a reads 

cj 01 = g'g dt , a/ 1 = g & , u lJ = Co lJ , (11) 
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The curvature £l ab of the connection uj ab reads 



n 2 V 



(12) 



2r 2r 
l_o3 . . 
O y = ^—E l E 3 . 



It is then straightforward to show that the entire content 
of ([HJ amounts to a single differential equation linear in 
g 2 — 1. Its solution is 



2 i 3ci 2 

9 - 1 = fe' 



c 2 A 
27^f 
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(13) 



Translating the c couplings into the more usual ones via 
(|4]) we indeed obtain the Boulware-Deser metric ([2]) with 
g 2 = Jbd- The integration constant C is of course re- 
lated to the mass of the solution. 

For future use let us note the following. Denote by 
gBD± the function g for the respective branch of the 
Boulware-Deser solution. One observes that 



2 i 2 

9bd+ "+" 9bd- 



3ci 

C2 



(14) 



This relation will be useful later on. 

The vacuum field equation 5C = for E is the match- 
ing condition 



i* s {- Cl e{(Q-u>)E 2 6E)+ 

+ -c 2 e((u> - + n - -(w - w) 2 }c5£)} = . 
2 3 



(15) 



if; denotes the pull-back of the form into E. 

Equation (fT5j) is derived under the hospitable assump- 
tion that the vielbein in spacetime is continuous across 
E modulo a Lorentz transformation. 

When E is time- or space-like, there is a unique geom- 
etry on E induced from the bulk. One may choose a viel- 
bein adapted to E [all but one its components are tan- 
gentially oriented] . Pulled back into E defines an induced 
vielbein, which is regarded as intrinsic to E. Solving the 
zero torsion condition one obtains a unique Levi-Civita 
spin connection for that vielbein. This is because pro- 
jecting r) a b into those tangential directions leaves us with 
an non-degenerate (invertible) tensor. The obtained con- 
nection coincides with the tangential components of lu. 
The induced fields from N and N may differ only by a 
Lorentz transformation. That is there is indeed a unique 
geometry on the hypersurface inherited from the bulk, 
and is regarded as intrinsic to it. Equation (|15[) . which in 
general involves a 'matter' energy tensor T a b on its r.h.s. 
according to ©, it is equivalent to the matching condi- 
tions first written down in 60] [6l| . In those works the 
matching conditions appeared in the usual kind of for- 
mulation (as set by the work of Israel 38] ) which involves 
explicitly the intrinsic geometry of E and the extrinsic 



curvatures of E w.r.t. the bulk regions N and N. The 
equivalence of that formulation to (fT5"|) was shown in [48| 
and re- visited in [57| . 

When E is null, the induced vielbein may also be re- 
garded as an intrinsic vielbein; it does span the space of 
tangential one-forms. The difference is that one direction 
is null and the projection of r\ a b into the tangential direc- 
tion is degenerate (non-invertible) . There is no unique 
Levi-Civita connection for the induced vielbein, thus no 
unique inherited geometry. 



41], see also [4£ 



Penrose 

classifies (three) different induced geometries which can 
be regarded as intrinsic to the hypersurface; these ge- 
ometries differ on the level the induced notions of paral- 
lel transport along E agree. These notions, presumably, 
depend on the derivatives of the metric. All this wealth 
of structure available in the null hypersurface is reflected 
in practice as follows. The normal vector is null therefore 
orthogonal to itself [and indeed tangential on the hyper- 
surface]. Thus one cannot construct projections along 
the tangential directions of E just by knowing its normal 
vector. That, combined with the non-uniqueness of the 
intrinsic geometry, implies that splitting the fields into 
components containing information intrinsic and extrin- 
sic to E is less straightforward in the null case. Formula- 
tions of the matching conditions involving explicitly the 
intrinsic geometry of the hypersurface necessarily have 
qualitative differences in the null and non-null case. Uni- 
fying descriptions of the two cases can be found [45] [47] 
but they rather emphasize the peculiarities of the null 
case. 

In our formulation of the matching conditions, leading 
to (|15[) . only bulk field information in the vicinity of E 
is involved. The intrinsic geometry information exists 
implicitly in the continuity conditions on the fields. That 
is, it exists in a minimal choice [so that to accommodate 
as many cases as cases as possible] of which field variables 
are continuous across E; the rest are discontinuous, as 
restricted of course by the field equations. Thus the 
causal character of E, time- or space-like, or null, enters 
the matching conditions only through our requirements 
on the fields in each specific problem. Formula ([15")) holds 
in all cases. 

In a most general setting, studied in Ref. 50], any 
chosen continuity conditions is such that the matching 
conditions are well defined, in the sense explained pre- 
viously. On this basis, the acceptable continuity con- 
ditions are certainly not unique [50]. In practice, one 
chooses acceptable continuity conditions which are ge- 
ometrically intelligible. As we have already mentioned, 
our choice here is a usual one: We have required that 
the veilbein in spacetime is continuous across E modulo 
a Lorentz transformation, which means that the metric 
tensor g — rj a t,E a (g> E b in spacetime is continuous. 

It is worth to mention that the choice of the vielbein 
E a and the connection uj a b as gravitational field vari- 
ables plays an important role. A choice of variables in- 
volving, for example, the metric tensor would necessar- 
ily involve its inverse. In any such case our formulation 
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of the problem in terms of globally defined discontinu- 
ous fields would become complicated and almost none 
of the tools used in Refs. |48T - l50l | could be applied effec- 
tively in order to derive equations such as (fT5j) and prove 
facts about them. Moreover, the one-forms E a and u a b 
are tensor valued in the Lorentz group. They can be 
written w.r.t. to arbitrary coordinates in each region 
in spacetime. The continuity condition on the veilbein 
across their common boundaries i.e. the hyper surf aces, 
is that the vielbein is continuous modulo Lorentz trans- 
formations, a local symmetry of the action (|5|). Indeed, 
in every specific problem, one can construct a Lorentz 
transformation between the veilbeins on the sides of a 
hypersurface. This transformation is then used in (jl5[) 
in order to express all fields in it w.r.t. the same basis; 
this completes the statement of the matching conditions. 
Therefore there is no need to explicitly use special coor- 
dinates chosen such that continuity is ensured. In this 
sense the formulation employed here parallels the work 
in Ref. ^7j> though the freedom to use any convenient 
coordinates on each side of the hypersurface is achieved 
quite differently. 

In all, one has in hand a formulation of the matching 
conditions which can be uniformly applied to the space- 
and time-like as well as to the null hypersurface which 
is the case of interest. Of course writing down the cor- 
rect equations of motion is different than understanding 
their content. In the null case the wealth of structure 
is revealed only after systematic general analysis. A full 
discussion of the null hypersurface matching conditions 
is far beyond the scope of the present work. In what 
follows we apply (1X5]) to show that the shock waves we 
advertised at the beginning of our work indeed exist. 



IV. SHOCK WAVES 

Continuity of the veilbein in spacetime M means the 
following. Each veilbein field, E and E, can be extended 
from the region it is defined, N and N respectively, 
across E and into a neighborhood of the other region. 
In the overlap E and E may differ by a Lorentz trans- 
formation which is a local symmetry of action (JSJ. Let 
us put this condition into formulas for the problem of 
interest. 

On the manifolds N and N the vielbein is given re- 
spectively by 

E°=gdt, El = ^J, E l =rE l (16) 



and 



Let S be a hypersurface along which 

dv dv — 

gdt = 0, —-gdt = 0, (18) 

9 9 

in N and N respectively. More precisely, the pull-back of 
the forms on the l.h.s. of these relations into S vanishes. 
E is null. 

A null hypersurface in spacetime dimension five is gen- 
erated by a three-parameter family of null geodesies, 
one through each point of the hypersurface. The null 
geodesies along E satisfy (fig)) , with parameters the 
points on the unit 3-sphere. Along these geodesies 
ds 2 = thus the induced metric is degenerate. The de- 
generate induced metric is another way, especially when 
working directly with the metric tensor, to recognize a 
null hypersurface. 

The field E a is assumed extendible across E and into 
the manifold N. Similarly the field E a is assumed ex- 
tendible across E and into the manifold N. Explicitly 
that means the following. There is a neighborhood of M 
around E where E a and E a co-exist. This is the over- 
lap of their support in spacetime. In that neighborhood 
define fields E a and E a such that 

E a is a rotated E a so that £ Q | S = £ a | s , (19) 
E a is a rotated E a so that E a \s — E a \x , 



By 'rotated' we of course mean Lorentz transformed 



E a = A a a E a , 



(20) 



everywhere in the overlap of E a and E a in M. The in- 
verse of the matrix (A a a ) is the matrix (A° a ). Therefore 
E a = A a a E a in that neighborhood. Completely analo- 
gous formulas hold for the unbarred fields as we extend 
E a into N, We shall not need explicitly these formulas. 

Continuity of the a — i components of the vielbein 
implies that r = f along E. That implies that 



^dr 



i^dr . 



(21) 



In view of ([T5jl . t + J drg~ 2 is a natural coordinate along 
E in the spacetime region TV. Similarly t + J dfg~ 2 in N. 
Therefore the forms dt+g~ 2 dr and dt + g~ 2 df may agree 
at E modulo suitable factors [null vectors have no natural 
normalization]. Indeed, setting g 2 dt + dr — g 2 dt + dr at 
E is consistent with ([2~Tj) . Translating that in vielbein 
components we write 



g{E° + E 1 )^ = g(E° + E 1 ) 



(22) 



It is then not hard to show that condition (|22|) implies 
a Lorentz transformation A a h which reads 



E u =gdt, 



E l 



dr 
9 



E l = fE l 



(17) 



Recall that E l is a veilbein on the unit round 3-sphere 
introduced in the previous section, formulas (j9]) and (fTQ|) . 



A a - = 




9 +9 



(23) 
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We suppress the trivial transformation of the angular 
components. 

This transformation had to be explicitly known in or- 
der to express all tensor valued forms in the same basis: 
the fields ui a r and Cl\ must be transformed to obtain the 

b b _ 

components of uj and O in the directions of the basis 
E a , which is the spacetime veilbein field appearing in 
equation (TT5|) . 

Under a Lorentz transformation A a a the spin connec- 
tion Co transforms as 

Cj\ = k\oj\k\ + k\dk\. (24) 

The transformation of its curvature tensor follows and 
reads 

n\ = k\n\k\. (25) 

These transformations hold in the overlap of the support 
of E and E in M. Under (f2"0"| and (f2"4")) the torsion tensor 
of ^transforms as dE a +ui a b E b = k a - a {dE a +ui\E l ). Be- 
ing indeed a tensor, its components in both bases vanish 
simulataneously. Therefore Hi a b , given by (|2"4"|). are indeed 
the components of the Levi-Civita connection Co in the 
basis which coincides with E a at E. Everything needed 
in order to apply (|15p has now been made explicit. 

The quantities of importance are the pull-backs of the 
'jumps' (Q— uj) a b into E. Recalling (TTTI) a straightforward 
calculation gives 

= 0, (26) 

-2 2 
•* f — \i i 9 9 friz 

Q 2 ~ Q 2 - ■ 
i e (w-cj) = ^— E . 

The jump of the angular components uj 1 ^ = u) % a vanishes 
trivially. It is then easy to calculate the pull-back into E 
of the jump squared ((a) — uj) 2 ) a b = (u> — u) a c {u — uj) b c , 
which vanishes identically. [Such a simplification, due 
to the high symmetry of the configuration, was vaguely 
conjectured in Ref. [62j|.] Using (TT21) the forms (fi + 0) a b 
are also straightforwardly calculated. 

Substituting everything in the matching condition 
(IT5]) one finally finds 

_ {sa+a L^£ } !^ x (2T) 

where J5~ := E° - E 1 . Also := + i? 1 . The 
vielbein (E ± ,E l ) is a basis adapted to E. 

E~ is the normal to E component of the vielbein in 
the sense that i^E~ = 0, by (fT5)) . The presence of 8E~ 
in (|27p is a reminiscent of the tangent nature of the null 
normal vector. Explicitly the variation reads 

i^SE- = \{xy + i'^E+ + X(x)~ j i* s E 3 , (28) 



for an arbitrary smooth field X(x) a b with support on E. 
The component of the matching condition (|27[) involving 
\(x)~ + it is not trivial. Therefore (|27p implies that 

(o , 2 -9 2 -g 2 \ g 2 - g 2 n , 0O ^ 

| 3ci+C2 _ j^^-O (29) 

for all r along E. The solution g 2 — g 2 amounts to no 
discontinuity at all. The non-trivial solution is 

Zci + c 2 2 ~ g2 r ~ g2 =0. (30) 

Recalling (|14|) we see that this equation can indeed be 
satisfied, and in fact g and g must belong to different 
branches of the same Boulware-Deser solution given in 
(J5J). Thus we obtain the branch changing shock wave so- 
lution we claimed to exist in the beginning of this paper. 

V. COMMENTS 

Mathematically, the spherical symmetric shock wave 
solutions of Lovelock gravity have been rather estab- 
lished. These shock waves are possible due to the multi- 
valued-ness of the spherically symmetric metrics of the 
theory. On the other hand, the multi-valued-ness itself 
is a symptom of a deeper pathology. 

The whole thing can be better understood considering 
a simpler system. An elementary such system is a point 
particle with Lagrangian L = fw +gv . v is the velocity 
of the particle and \i and g are non-zero constants. The 
canonical momentum reads p = 2^lv + Agv 3 . Varying the 
action J dtL one obtains the Euler-Lagrange equation, 
p = 0, the dot denoting time derivative. That is, p 
is constant in time. Let across some instant of time 
momentum change from p ln to p ou t and velocity from Vi n 
to v ou t- The equation of motion requires that pi n — p ou t- 
This translates to 

(fin - v out )(2fi + ig(vf n + v in v out + w„ ut )) = . (31) 

This is a matching condition. Along with the solution 
Vin — Vout = i.e. velocity is constant, one obtains a solu- 
tion in which velocity jumps according to the condition 
vf n + WmWout + v 2 out = -[i I (2g). The jumps may happen 
at arbitrary times. We learn that a source-free motion is 
not necessarily a uniform motion in this system. In fact, 
p = constant may correspond to an infinity of different 
piece- wise uniform motions. 

From the Hamiltonian point of view, the problem 
arises in the following form. The Hamiltonian is a 
function of the canonical momentum p. Solving p — 
2(xv + Agv 3 for v the answer is multi- valued. Each one 
of the multiple solutions v(p) defines a different 'canon- 
ical branch' of the theory, as we may call it, with a 
different Hamiltonian H = pv(p) — L. We have seen 
that the source-free motion p — constant is neither sim- 
ple nor unique in this theory, and in general involves 
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sudden changes of velocity according to the equation 
vj n + Vj n v out + Vo Ut = —fi/(2g). These jumps are a sudden 
change of canonical branch as vi n and v ou t are values of 
the velocity given by different v(p). Thus the Hamilto- 
nian by which the system evolves may change abruptly 
without a cause. A byproduct of this fact is that energy 
is not conserved in the free motion of this system. More- 
over, one cannot uniquely relate the space of solutions 
to the space of the initial data i.e. there is no notion of 
classical phase space for this theory. 

All that translates as it is in Lovelock gravity. The 
role of velocity is played by the extrinsic curvature which 
jumps across a hypersurface in vacuum. The Lagrangian 
is not linear in the curvature i.e. not quadratic in the ex- 
trinsic curvature. This leads, in the quise of many com- 
ponents and complicated formulas, to the same problems 
encountered in the simple point particle system. The 
multi-valued-ness of the Lovelock gravity Hamiltonian 
and some of its implications have been emphasized in 
Ref. [63 |. In Ref. 37 1 it was shown by explicit examples 
that appropriate jumps of the extrinsic curvature across 
non-null hypersurfaces are allowed in vacuum, i.e. the 
canonical momentum does not change across the hyper- 
surface nonetheless the theory changes canonical branch. 
Mathematically the results arise through factorized ex- 
pressions analogous to equation (|3"Tj) : There is a continu- 
ity imposing factor which we require not to vanish, and 
another factor originating from the higher order terms 
which leads to a soluble equation. That was realized 
here in the form of equation (|29p leading to equation 
(|3U|) and our shock wave solutions. The factorization 
can be seen already at the level of equation (fT5|) which 
is a complicate analogue of (|3l"j) . In fact factoring out 
Q — lu in (|15l) one obtains a rank three tensor whose van- 
ishing is the general condition that must be satisfied by 
the fields in a change of canonical branch. The multi- 
ple canonical branches are due to the non-linearity of 
the Lagrangian in curvature, an outcome of which is the 
multi-valued-ness of solutions. 

The shock waves presented here are spherically 
symmetric and exist in dimension higher than four. 
It takes the peculiar dynamics of Lovelock gravity for 



them to arise. Now it is worth to mention that, even 
in Einstein gravity, once in higher dimensions spherical 
symmetry is not a necessary ingredient of Birkhoff's 
type of theorems. The reason why is that Israel's black 
hole uniqueness theorem [64j [65] , which says that static 
and asymptotically flat black holes are necessarily spher- 
ically symmetric, is not as strong in higher dimensions 
as it is in four: If one drops asymptotic flatness then the 
geometry of the 'angular' manifold [with the topology of 
the sphere] is not necessarily that of the round sphere, 
see the discussion in [66j[67|. Einstein field equations 
restrict the 'angular' manifold geometry only that much. 
Now all that 'room to spare' might make possible even 
time- dep endent solutions as it is numerically argued in 
Ref. [68j using an angular manifold with the topology 
of the three-sphere but not the geometry of the round 
sphere. In fact all that freedom is controlled once one 
switches on the Lovelock gravity terms available in those 
higher dimensions, as emphasized in the Refs. [69j][36j]. 
In this work we have seen that in Lovelock gravity C° 
piecewise smooth time-dependent solutions exist in the 
form of shock waves even under spherical symmetry. 
The derivation goes through as it is with minor modifi- 
cations in any other case of interest [one simply uses an 
'angular' vielbein E l with the appropriate properties] 
e.g. the Boulware-Deser-Cai metrics [70] or the metrics 
of Ref. (3(| solutions of theory (fTJ). Formulas become 
only more cumbersome, not essentially different, when 
including cubic or higher Lovelock terms in the appro- 
priate dimension. The analogues of the shock waves 
presented here arise as long as the smooth solutions 
are multi-valued through elementary identities such as 
relation (fT4l). 
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